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Abstract 

We present a new high-order coupled cluster method (CCM) formalism for the ground states 
of lattice quantum spin systems for general spin quantum number, s. This new "general-s" 
formalism is found to be highly suitable for a computational implementation, and the technical 
details of this implementation are given. To illustrate our new formalism we perform high-order 
CCM calculations for the one-dimensional spin-half and spin-one antiferromagnetic XXZ models 
and for the one-dimensional spin-half /spin-one ferrimagnetic XXZ model. The results for the 
ground-state properties of the isotropic points of these systems are seen to be in excellent 
quantitative agreement with exact results for the special case of the spin-half antiferromagnet 
and results of density matrix renormalisation group (DMRG) calculations for the other systems. 
Extrapolated CCM results for the sublattice magnetisation of the spin-half antiferromagnet 
closely follow the exact Bethe Ansatz solution, which contains an infinite-order phase transition 
at A = 1. By contrast, extrapolated CCM results for the sublattice magnetisation of the 
spin-one antiferromagnet using this same scheme are seen to go to zero at A « 1.2, which 
is in excellent agreement with the value for the onset of the Haldane phase for this model. 
Results for sublattice magnetisations of the ferrimagnet for both the spin-half and spin-one spins 
are non-zero and finite across a wide range of A, up to and including the Heisenberg point at A = 1. 

Running Head: CCM Calculations of Quantum Magnets with General Spin 
Keywords: CCM, Quantum Magnets, Phase Transitions. 
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I. INTRODUCTION 



In this article we consider lattice quantum spin systems at zero temperature. Since the 
inception of this subject via the introduction of the Heisenberg model it has undergone a 
period of steady development. Although the exact Bethe Ansatz solution of the Heisenberg 
model for the spin-half linear chain 

mm 

followed quickly after its initial introduction, 
interest in lattice quantum spin systems has not faltered ever since. In particular, exact 
solutions, such as those typified by the Bethe Ansatz, are generally restricted to unfrustrated 
systems of low-dimensional lattices and to low spin quantum number. Thus, much effort 
has been expended over the last fifty or so years in trying to understand the properties of 
more general lattice quantum spin systems at zero temperature via approximate methods. 

A recent and exciting topic in this field has been the development of the density matrix 
renormalisation group (DMRG) 

mm 

method. This technique allows one to per- 
form highly accurate calculations for both frustrated and unfrustrated systems of general 
spin quantum number. As long as one is able to partition a one-dimensional quantum sys- 
tem into "system" and "environment" subsystems then one ought in principle to be able 
to apply this powerful method of modern-day quantum theory. In particular, DMRG cal- 
culations have been used to demonstrate conclusively that the spin-one Heisenberg linear 
chain antiferromagnet contains an excitation gap 10|. This is in stark contrast to both 
its classical behaviour and the behaviour of its quantum spin-half counterpart. However, a 
major restriction on the DMRG method is that it has, so far, only been conclusively ap- 
plied to one-dimensional or quasi-one- dimensional systems, although isolated cases of highly 
successful DMRG calculations for various two-dimensional lattices have been performed. 

By con.a,.. ,ua,.u. Mo,.e CaHo ,QMC) ca.cu.a.o„. Q fl U 13 have bee„ a.- 
ready applied with much success to spin systems of spatial dimensionality both equal to and 
greater than one. In particular, QMC calculations for the zero-temperature properties of 
the spin-half, square-lattice lo, ll^ ll^ Heisenberg model present a very accurate and 
valuable benchmark against which other approximate methods may test themselves. How- 
ever, although QMC is not as limited by spatial dimensionality as the DMRG method, it is 
limited in its range of applicability by the existence of the well-known "minus-sign problem." 
For lattice quantum spin systems the minus sign problem is, in turn, often a consequence of 
frustration. We note however that for non-frustrated systems one can often, but not always. 



3 



determine a "sign rule" 
QMC. 

Other approximate methods that have been apphed to 
zero temperature) include spin-wave theory (SWT) 



29| which effectively circumvents the minus-sign problem of 
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59|. For SWT one maps the spin system onto an exact 
(or effective) bosonic system (for example, via the Holstein-Primakov transformation) which 
may then be solved. The advantage of linear SWT is that it is simple to apply and an exact 
solution generally exists via the Bogoliubov transformation. Higher-order SWT becomes 
rapidly more complicated and may only be solved perturbatively for these higher orders. 
Also, it is difficult to increase the level of approximation of this method in a structured and 
well-defined manner. By contrast, for the QMC method one may systematically increase 
the lattice size and still obtain an "exact" result to within a well-defined statistical error. 
Exact diagonalisations of finite-sized systems [2^ 24 1 have also been applied to quantum 



spin models with much success. However, such calculations rapidly become limited as the 
size of the lattice is increased by the amount of the computational resources available. &\ 
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contrast, this problem is not so evident for the exact cumulant series expansion 
calculations. Indeed, such calculations have provided very accurate results for the ground- 
and excited-state properties of various spin systems. However, in order to determine such 
expectation values a resummation of the otherwise rapidly divergent perturbation series 
must be performed using Fade approximants or related techniques. Finally, a technique of 
quantum many-body theory (QMBT) called the coup 
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so been applied 

to both unfrustrated and frustrated lattice quantum spin 
systems with considerable success, and it is this method upon which we shall concentrate 
henceforth. 

The first application [3l of the CCM to the lattice quantum spin systems was to the spin- 
half Heisenberg model and related models on various bipartite lattices using several localised 
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4l| to include anisotropy 



approximation schemes. This work was quickly generalised 
within the Hamiltonian via the spin-half XXZ model. In particular, critical points as a 
function of the anisotropy were discovered in an approximation scheme denoted as SUB2 
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(in which all possible two-spin correlations are retained). It was furthermore found at these 
critical points that the CCM excitation spectra became "soft" and that the behaviour with 
respect to spatial separation of the spin-spin correlation functions changes from exponential 
to algebraic decay. Both of these features were clearly strong indicators that these critical 
points were reflections of quantum phase transitions in the real system, which in turn were 
driven by the anisotropy parameter within the Hamiltonian. (An interesting application 
of the extended coupled cluster method (ECCM) to this model was also performed, and 
the interested reader is referred to Ref. [5J|.) Frustrated quantum spin systems, such as 
the J1-J2 model and the triangular lattice antiferromagnet have 

also been studied using the CCM with equal success. Indeed, the CCM has been utilised 
not only to provide accurate values of the ground-state properties of these systems and 
very accurate predictions of their phase diagrams, but also to simulate the nodal surfaces 
I52L l53| of these models. We note that such approximately determined nodal surfaces are 
potentially of very great use to QMC calculations in which the "sign problem" is p resent. 
Other such applications of note of the CCM have been to the biquadratic model and to 
the one-dimensional spin-half J1-J2 model using a "dimerised" model state l44l . 
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Furthermore, recent high-order CCM calculations j4 
57, 5^ for spin-half systems using localised approximation schemes have been shown 
to be very effective in simulating the ground- and excited-state properties of quantum spin 
systems. The strength of this approach is that one can systematically and rigorously increase 
the approximation level as a function of some parameter which increases the size of the 
"locale" over which one retains all multi-spin correlations. The approximation schemes 
become exact in the limit of this parameter going to infinity, although this is generally 



impossible to achieve in a practical application. One can however extrapolate [55| a series 
of such results to this limit, and thus obtain very accurate approximate results for the 
exact properties of even highly frustrated systems. Indeed, the results produced for both 



unfrustrated systems such as the spin-half XXZ |46 . 



3, HQ, 



Ising models 



and frustrated systems such as the spin-half J1-J2 model 



and transverse 
Q|, the spin- 



half triangular [5l| and kagome '58] lattice antiferrornagnets, and the spin-half square-lattice 
Heisenberg antiferromagnet with two types of bonds 53|, have all been seen to be in excellent 
agreement with the results of the best of other approximate methods, where such results 
exist. A recent calculation has been performed for the spin-one square-lattice XXZ model 
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59| which utihsed the new "general-s" formahsm presented in this article, and the results 
for the ground-state properties of this model were found to be in superb agreement with 
previous SWT and cumulant series expansion results. 

We note that although spin problems are conceptually simple, they often demonstrate rich 
and unusual phase diagrams which are the direct result of the strong influence of quantum 
fluctuations in these strongly correlated systems. Indeed, these phase diagrams can become 
richer for systems with spin quantum, s, greater than 1/2, although we note that the limit 
s — > oo coincides with the pure classical model. For example, the spin-one biquadratic 
model contains antiferromagnetic, dimerised, and trimerised phases, and a "Haldane" 
phase as a function of the ratio of the bond strengths of the linear and biquadratic terms. 
Thus, lattice quantum spin problems open a wide window on to the rapidly developing field 
of quantum phase transitions. This subject becomes even more interesting when one also 
considers the non-zero temperature behaviour of such systems, and the respective roles of 
quantum and thermal fluctuations in driving phase transitions. For the present, however, 
we restrict ourselves to the zero-temperature case. 



II. THE COUPLED CLUSTER METHOD (CCM) 
A. The Ground- State Formalism 

The exact ket and bra ground-state energy eigenvectors, and of a many-body 
system described by a Hamiltonian H, 

H\^)=E,\^); {m = E,{^\, (1) 
are parametrised within the single-reference CCM as follows: 

|vl/)=e^|$) ; S = Y.SiCi, 
(§| = ($|5e-^ ; S = l + J2SjCT. (2) 

The single model or reference state |$) is normalised (($|$) = 1), and is required to have 
the property of being a cyclic vector with respect to two well-defined Abelian subalgebras 
of multi-configurational creation operators {C/} and their Hermitian-adjoint destruction 
counterparts {Cj = (C/)"''}. Thus, |$) plays the role of a vacuum state with respect to a 
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suitable set of (mutually commuting) many-body creation operators {C/}, 

crm = o , j^o, (3) 

with Cq = 1, the identity operator. These operators are complete in the many-body Hilbert 
(or Fock) space, 

We note that although the manifest hermiticity, ((^^|^ = |\E')/(\1/|\I')), is lost in these 
parametrisations, the intermediate normalisation condition (^E'l^E') = (*&|^) = ('^1'^') = 1 
is explicitly imposed. The correlation coefficients {Si,Si} are regarded as being indepen- 
dent variables, even though formally we have the relation 

(<l>|e^^e^ 

{^\s= ' t ^1^, ■ (5) 

The full set {Sj, Sj} thus provides a complete description of the ground state. For instance, 
an arbitrary operator A will have a ground-state expectation value given as 

A= {^\A\<il) = {^ISe-^Ae'^l^) = A(^{Sj,Sj}) . (6) 

We note that the exponentiated form of the ground-state CCM parametrisation of Eq. (j21) 
ensures the correct counting of the independent and excited correlated many-body clusters 
with respect to |$) which are present in the exact ground state |\1'). It also ensures the 
exact incorporation of the Goldstone linked-cluster theorem, which itself guarantees the 
size-extensivity of all relevant extensive physical quantities j37|. 

The determination of the correlation coefficients {Sj,Sj} is achieved by taking appropri- 
ate projections onto the ground-state Schrodinger equations of Eq. ((H). Equivalently, they 
may be determined variationally, 

SH/6S1 = ^ ($|C7e-^i/e^|$) = 0, V/ ^ ; (7) 
SH/SSi = mSe-'^[H,C+y\<^) =0, V/ 7^ . (8) 

Equation (|7j) also shows that the ground-state energy at the stationary point has the simple 
form 

E, = E,i{Sj}) = {^\e-'He'\^) . (9) 
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We note that Eq. ((Tj) represents a coupled set of nonlinear multinomial equations for the 
c-number correlation coefficients {Sj}. The nested commutator expansion of the similarity- 
transformed Hamiltonian, 

H ^ e-'He"^ = H+[H,S] + ^[[H, S], S] + ■ ■ ■ , (10) 

together with the fact that all of the individual components of 5* in the sum in Eq. (j21) 
commute with one another, imply that each element of S in Eq. is linked directly to the 
Hamiltonian in each of the terms in Eq. (fTIUl . Thus, each of the coupled equations (jZj) is of 
linked-cluster type. Furthermore, each of these equations is of finite length when expanded, 
since the otherwise infinite series of Eq. fllO|) will always terminate at a finite order, provided 
(as is usually the case) only that each term in the second-quantised form of the Hamiltonian 
H contains a finite number of single-body destruction operators, defined with respect to 
the reference (vacuum) state 1$). Therefore, the CCM parametrisation naturally leads to a 
workable scheme which can be efficiently implemented computationally. 

The CCM formalism is exact in the limit of inclusion of all possible multi-spin cluster 
correlations within S and S, although in any real application this is usually impossible to 
achieve. It is therefore necessary to utilise various approximation schemes within 5* and S. 
The three most commonly employed schemes have been: (1) the SUBn scheme, in which all 
correlations involving only n or fewer spins are retained, but no further restriction is made 
concerning their spatial separation on the lattice; (2) the SUBn-m sub-approximation, in 
which all SUBn correlations spanning a range of no more than m contiguous lattice sites 
are retained; and (3) the localised LSUBm scheme, in which all multi-spin correlations 
over distinct locales on the lattice defined by m or fewer contiguous sites are retained. 
We also make the specific restriction that the creation operators {C/} in S preserve any 
additional symmetries of the Hamiltonian. Thus, the approximate CCM ground-state wave 
function is constrained to lie in the appropriate subspace defined by the additional quantum 
numbers corresponding to these additional symmetries. For the XXZ model illustrated later 
the additional symmetry is provided by the total ^-component of spin sf. = J^i^i, which 
commutes with the Hamiltonian. The ground state lies in the sector = 0. We denote as 
distinct configurations those in such appropriately defined subspace which are inequivalent 
under the point- and space-group symmetries of both the lattice and the Hamiltonian. The 
number of such distinct (or fundamental) configurations for the ground state at a given level 
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of approximation is labelled by Np. 



B. The High-Order Formalism For General Quantum Spin Number 

In order to determine the CCM ground-state ket configurations we fundamentally need 
to pattern-match the configurations in the set {C/} to the spin-raising operators contained 
in H\^). We note that for small values of the truncation indices {m,n} mentioned in the 
previous section (and thus for low orders of approximation) this may readily be performed 
analytically. However, for higher orders of approximation we must use computational meth- 



ods (see for example also Refs. 46|, |47|, |5l|, |55|) in order to do this. For the cases of interest 



here, we begin by defining a set of local spin axes in which all of the spins in the chosen 
model state |$) point along the respective negative z-axes, namely 

N 

1$) = I l)i ; in the local quantization axes, (11) 

1=1 

where | |)j = |s,— s)j. This is achieved by an appropriate set of local rotations. Since 
such rotations are canonical transformations, the underlying spin algebra is preserved, and 
the energy spectrum of the transformed Hamiltonian (i.e., written in the rotated local spin 
coordinate scheme) is unchanged. 

The next step in the computational implementation of the CCM for lattice quantum spin 
systems of general spin quantum number, s, is to define a suitable set of multi-spin creation 
and destruction operators with respect to this model state. We thus define the CCM ket- 
state correlation operator S in terms of sums of products of single spin-raising operators, 
s'l = sf. + isl, (again with respect to their local spin axes), such that 

N N 

S = J2S.,sl + Y.S^„^,st,sZ + ■■■ . (12) 

The coefficients iSj^, Si^^^i^, and so on, now represent the spin-correlation coefficients specified 
by the sets of site indices, {^i}, {ii, ^2} and so on, on the regular lattices under consideration. 
We note that these indices run over all lattice sites, and that different indices may thus 
indicate the same lattice site. For the case of general spin quantum number s we note that 
we have a maximum number of spin-raising operators at any specific site / which is 2si, 
where s; is the spin quantum number of the spin situated at site /. For the spin-half case, we 
are thus limited to only one spin-raising operator per lattice site as required, and we note 



ihat in this manner we build in a previous spin-half high-order CCM formalism (see Refs. 



46 



47 



51 



55|) into the new high-order formalism for general-s directly from the outset. In 



order to simplify the high-order CCM formalism, it is also found to be useful to define the 
following operators: 



El E 



Fk 

Mk^rn,n = El>2 X]j4, ~ ^ 2) <Sk,m,n,i4,---,ii ^ 

Nk,m,n,p = j:i>3 Ei^-M i{i-m-m-'i)Sk 



■m,n,p,i5,---,ti '^45 



(13) 



Hence, we determine the similarity-transformed expressions of the single-spin operators 
3°" ; a = {+,—,z}, where = {s^Y = — is^, by using Eq. (jH]) and the usual spin 



commutation relations, such that 



Si. e" 



e ^sle^ = si 



Su e" 



St + Fksj 



'^Fusl 



Gkk^k 



(14) 



For the specific case of s = 1/2 at site k we note that Gkk = because "double occupancy" of 
the lattice site k is prohibited in this case. In order to determine the similarity transformed 
version of a given Hamiltonian, we also need to know the commutation relations of the 
operators defined in Eq. (fTSj) with the single-spin operators s" ; a = {+, —,z}, and these 
are stated in the Appendix. 

We now define the set of CCM destruction operators {Cf} (of / number of spin-lowering 
operators), as follows 

^~ ~ ~ s. 



where the indices ji,j2, ■ ■ -yji represent any given lattice site. We choose only one of the 
NbQ'^^^i symmetry-equivalent configurations to pattern-match with the terms within H in 
order to determine the CCM ket-state equations of Eq. (|7j). Note that A''^ is the number of 
Bravais lattice sites and that, for a given cluster /, vi is a symmetry factor dependent purely 
on the point-group symmetries (and not the translational symmetries) of the crystallographic 
lattice. 

The process of the enumeration of all possible fundamental clusters and the process of 
"pattern-matching" are both ideally suited to an efficient computational implementation. 
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and a full description of these processes is also given in the Appendix. Furthermore, it is 
a simple matter to determine and solve the CCM bra-state equations, once the ket-state 
equations have been obtained and solved, as described in the Appendix, where we also 
explain the technicalities involved in obtaining ground-state expectation values. 

In the remainder of this paper we now demonstrate the new formalism for various models 
of interest. 



III. THE ONE-DIMENSIONAL HEISENBERG MODEL 

The ID XXZ model is defined by the following Hamiltonian, 

E = j:{s^s^^, + sysl, + AstslA , (16) 
i=i ^ ^ 

where the index i in Eq. runs over all lattice sites on the linear chain. We are 

interested specifically in the infinite chain, A^ ^ oo. The system is by now extremely 
well understood for the spin-half case via the existence of the Bethe Ansatz 
exact solution for this case. The exact ground-state energy for the infinite chain at A = 1 
was shown to be given by the expression, Eg/N=l/A + ln(l/2) —0.44314718). Indeed, 
the long-range Neel ordering inherent in the solution in the limit A ^ oo is completely 
destroyed by quantum fluctuations at the phase transition point at A = 1. By contrast, no 
exact solution for the spin-one Heisenberg model on the infinite linear chain has as yet been 
found, although extremely accurate DMRG calculations have been performed [8] for this 
model giving a value for the ground-state energy of £'g/A^=— 1.401484038971(4) at A = 1. 
Furthermore, these calculations conclusively showed that there is an excitation energy gap 
of magnitude 0.41050(2) in this system, which was first postulated by Haldane This 
is in stark contrast to the spin-half model which is gapless, and we note that conventional 
spin-wave theory predicts that the excitation energy of the Heisenberg model on the infinite 
linear chain for both the spin-half and spin-one cases is gapless. For the spin-one anisotropic 
Heisenberg model, the isotropic Heisenberg point is in the Haldane phase, in which the 
amount of long-range Neel-ordering is zero and the excitation spectrum is gapped. The 

HLse transition from the Neel-like phase occurs for a value of anisotropy A = 1.167 ± 0.007 
for the spin-one XXZ model. 
Interest has also been expressed recently in one-dimensional ferrimagnetic Heisenberg 
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systems, for which case the Lieb-Mattis theorem states that the system may demonstrate a 
macroscopic lattice (and not sublattice) magnetisation. For example, we note that for the 
spin-half /spin-one ferrimagnets in which alternating spins on the chain have, respectively, 
spin quantum numbers of 1/2 and 1, recent DMRG calculations jol] predict that the ground- 
state energy of the Heisenberg model (A = 1) is given by Eg/N ^ —0.72704 and the 
sublattice magnetisation on the spin-one sublattice is given by 0.79248 and the sublattice 
magnetisation (multiplied by a factor of 2 here - see below) on the spin-half sublattice is 
given by 0.58496. The phase transition is also believed to be at (or near to) the Heisenberg 
point (A = 1). 

The model state that we shall use in our COM calculations here is the one-dimensional 
Neel state. We perform a rotation of the local spin axes of 180° about the y-axis of the spins 
on one sublattice only, such that the model state now appears mathematically to consist 
entirely of spins which points along the negative 2;-axis, as in Eq. ()11|) . We note however 
that this unitary transformation does not affect the eigenvalue spectrum of this problem. 
The Hamiltonian may now be written in terms of these local spin axes, as 

H = - e{ A^.^^^+i + 1 {4sU + s-^r+i) } , (17) 

where = s^±is^. For the spin-half and spin-one antiferromagnets, we define the sublattice 
magnetisation in terms of these local spin-axes as, 

1 ^ ~ 

M = --T^T.i'^lsm , (18) 

1=1 

where the index i runs over all N lattice sites. This reflects the fact that for the spin- 
half and spin-one Heisenberg models we have a unit cell which contains one single lattice 
site only. By contrast, for the spin-half /spin-one ferrimagnet we have a unit cell which 
contains two nearest-neighbour spins of quantum spin numbers 1/2 and 1. We thus define 
the sublattice magnetisations Mi, on the A (spin-half) sublattice sites, and M2, on the B 
(spin-one) sublattice sites, separately, in the local spin axes by 

2 ~ 

Ml = -TrE(^l4l*) and (19) 

1 ii 
1 ^2 ^ 

where ii runs over all A^^i = N/2 spin-half sites and ^2 runs over all N2 = N/2 spin-one sites. 
Furthermore, we note that the Lieb-Mattis theorem states that \Mi/2 — M2I = 1/2 (after 
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rotation of the local spin axes) thus allowing an overall macroscopic magnetic moment. The 
technicalities of determining M, Mi, and M2 are explained in the Appendix. 

We may now perform the similarity transform of the Hamiltonian in order to write that 
Hamiltonian, in terms of the operators F^, Gk^m, Mk,m,ni and Nk^jn,n,p given above. By 
making use of Eqs. (HH) and (jMl) we find that H\^) = e'^He'^l^) = {Hi + H2 + ^3)!'^'), 
where 

+4FjFj+iG'j,i+i + Gi^iGi+i^i+i + Gi^iFl_^_^ + F^Gi+i^i+i + F^FIj^^ | s^^s/^i , (21) 
H2 = - Y.[UF., + M,,,^i + 2F,G,,+i + G,,F.,^i + FfF,+Xtsl^^ 

(22) 

^3 = -Ei^ + SGM+i + Si^F.+ij . (23) 



Again, we note that repeated indices in the operators Gk,m, Mk,m,n, and Nk,m,n,p operators are 



prohibited for the spin-half case. Hence, the general spin quantum number forma 
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ism reduces 
51, y. A 



to the previous formalism determined for the spin-half case in this limit 
general feature of CCM calculations for spin systems is that we often solve the CCM at 
different values of some tunable parameter within the Hamiltonian. For the case of the 
Hamiltonian of Eq. p7|l this is the anisotropy parameter. A, and we note that we "track" 
the solution from the trivial limit A ^ 00 to lower values of A. Another feature of CCM 
calculations is that often the real solution to the CCM equations is seen to terminate at 
some critical value of this parameter (e.g., denoted Ac here), and this behaviour is taken to 
be a signal of a phase transition in the real system. This illustrates a strong advantage of 
this method: namely, that the CCM is able to make predictions for the positions of quantum 
phase transition points from within a fully ah initio framework. 



IV. RESULTS 



In this article we utilise two approximation schemes in order to determine CCM expec- 
tation values for the ground-state energy and the sublattice magnetisation. The latter gives 
a measure of the amount of sublattice ordering. The two schemes are, namely, the LSUBm 
and SUBm-m approximation schemes, both of which include correlations in a systematic 
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and structured manner and provide exact results in the asymptotic limit m oo. We note 
that the LSUBm and SUBm-m schemes are equivalent for the ferrimagnet, if the truncation 
index m is an even number. This is due to the restriction that = J2i = in the ground 
state thus ruling out any clusters in the LSUBm approximation which contain more than m 
spin flips. This is noi the case for the spin-one antiferromagnet, in which case the restriction 
= J2i sf = may be satisfied for various clusters in LSUBm which do contain more than 
m spin flips. 

However, we note that it is usually impossible in a practical application of this method 
to exactly solve for the LSUBm and SUBm-m schemes in the limit m — ^ oo, and so we 
extrapolate the "raw" results as a function of m in this limit. Although no rigorous ex- 
trapolation scaling laws exist, we may attempt to extrapolate these results "heuristically." 
Thus, for example, extrapolation of these results in the limit m oo at each value of A 
independently may be achieved by assuming a "power-law" dependence, hereforth denoted 
as the "extrapolated CCM" results, given by 

yi = a + bx'^ , (24) 

where t/i is a CCM expectation value at a given level of LSUBm or SUBm-m approximation 
and Xi = 1/m, and there are p data points (i.e., for the p different values of m for which 
calculations have been performed). The best fit of this data to the parametric form given in 
the equation above is found and the extrapolated result is thus given by a. The interested 
reader will find a comprehensive explanation of the extrapolation process of CCM LSUBm 
expectation values for the spin-half XXZ model for a variety of lattices in Refs. [ssl . Isoj l . 

Results for the XXZ model are presented in Figs. CHSl and results for the Heisenberg 
model (A = 1) on the linear chain are given in Tables HI Hllll We may see from Figs. nand|21 
for the spin-half and spin-one XXZ models that our results for the ground-state energy are 
highly converged, and analogous behaviour is seen for the ground-state energies of the spin- 
half /spin-one ferrimagnet. Indeed, we may see from Tables HHITTl that our CCM results for 
the Heisenberg model (A = 1) converge rapidly in each case, and that they agree well with 
their respective exact or DMRG results. We note that the extrapolated CCM results for the 
spin-one Heisenberg model seem to agree with DMRG results to only two decimal places, 
whereas the extrapolated CCM results for the spin-half Heisenberg model appear to agree 
with exact Bethe Ansatz results to about five decimal places. The apparent discrepancy is 
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understood by noting that we are well into the Haldane phase at A = 1 for the spin-one 
model and the fact that we obtain results at this point at all is a testament to the power of 
the CCM. We note however that other model states might be employed at the Heisenberg 
point for the spin-one model in order to obtain even more accurate results for this model. 
Inspection of Table |llll f for A = 1) indicate that our "raw" (unextrapolated) CCM results for 
the ground-state energy of the spin-half/ spin-one ferrimagnet appear already to be converged 
to at least five decimal places. We expect that the CCM results are even more accurate for 
A > 1 because quantum fluctuations are (comparatively) weaker in this regime than at the 
isotropic Heisenberg point at A = 1. We furthermore expect that our extrapolated CCM 
value at A = 1 given in Table lllll represents an even better result than the raw SUBm-m 
results. Indeed, we note that all of our CCM results are in excellent agreement with those 
result of the DMRG method at this point. 

Figs. El and |3] present results for the sublattice magnetisation, M, for the spin-half and 
spin-one XXZ models in the Neel-ordered regime. The extrapolated CCM results for the 
spin-half XXZ model are clearly in excellent agreement with exact results and we note in 
particular that the sublattice magnetisation goes to zero at A = 1.07. This is in good 
agreement with the exact behaviour of the model and we note that there is an infinite-order 
phase transition at A = 1 to disordered phase for A < 1 in the 'real' system. Similarly, the 
extrapolated CCM results for both the SUBm-m and LSUBm approximation schemes for the 
spin-one XXZ model go to zero at A ^ 1.19 and A ^ 1.18, respectively. This is in excellent 
agreement with the position for the onset of the Haldane phase of A = 1.167±0.007 We 
furthermore note that even greater accuracy for the CCM extrapolated results for both the 
spin-half and spin-one antiferromagnets is expected with higher levels of approximation. We 
also note that these results are a clear indication that the CCM is capturing the qualitative 
difference in the behaviour of the sublattice magnetisations between the spin-half and spin- 
one antiferromagnetic XXZ models. 

Furthermore, an overall magnetic moment is possible for the case of the spin-half/spin- 
one ferrimagnet and we may see from Fig. and Table IIIII that the CCM results for the 
amount of sublattice ordering on both sublattices, namely. Mi and M2, are extremely well 
converged over the whole of the Neel-like regime. Indeed, we note that even though the 
values for Mi and M2 are determined separately using the CCM, they still obey relationship 
\Mi/2 — M2I = 1/2 specified by the Lieb-Mattis theorem. It is furthermore not a priori 
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evident that this ought to be the case as we have not artificially constrained this to be true 
in our CCM calculation at any point. It is thus a reflection that the CCM is detecting 
the underlying nature of the ground state from within an ab initio framework. We also 
see from Table UTTl that the CCM results agree with DMRG results to (at least) about five 
decimal places for the isotropic model. Furthermore, we believe that the extrapolated CCM 
results for the sublattice magnetisations at the Heisenberg point at A = 1 present even 
more accurate estimates of these quantities. We note that the results for the sublattice 
magnetisations of the spin-half /spin-one ferrimagnet remain non-zero and finite over the 
whole of the regime A > 1, although it is possible that a phase transition to a regime in 
which spins lie in the xy-plane (and for which Mi = M2 = 0) occurs for A ^ 1 (for example, 
see Refs. (We note that we would probably employ a model state in which spin 

lie in the xy-plane for the CCM in the region A < 1.) Our highly converged CCM results 
thus indicate that if the phase transition does indeed occur near to A 1 then this might 
involve a stepwise change in the sublattice magnetisations from a non-zero (and finite) value 
to zero at this point. If this conjecture were true then this behaviour would furthermore be 
consistent with a first-order phase transition at this point. 

Finally, we note that the CCM critical points, A^, are observed for the spin-one antifer- 
romagnetic XXZ model (illustrated in Table ITT]) at which the real (i.e., physical) solution 
to the CCM equations breaks down at given level of LSUBm or SUBm-m approximation 
level. This behaviour is an indication that a phase transition occurs in the "real" system. 
We note that an alternative and independent indication of the position of the phase tran- 
sition is afforded by the value of Ac at which the sublattice magnetisation goes to zero 
(discussed above). No equivalent results for the critical points of the ferrimagnetic XXZ 
model were obtained, although CCM calculations for this model show characteristic peaks 
in the derivatives of the ground-state energy for A < 1. This behaviour might possibly 
again be interpreted as a signature of a phase transition. It is also quite possible however 
that this typical breakdown of the real solution of the CCM equations might however occur 
at higher orders of approximation j48| . 
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V. CONCLUSIONS 



In this article a high-order CCM ground-state formahsm for lattices comprised of spins 
with general spin quantum number s was presented. We have shown that this new formalism 
is highly suitable for computational implementation for localised approximation schemes 
(namely, the LSUBm and the SUBm-m approximation schemes). In order to demonstrate 
the effectiveness of the formalism, successful applications of the general-s formalism were 
made to the one- dimensional XXZ model for the spin-half and spin-one antiferromagnetic 
cases and the spin-half/spin-one ferrimagnetic case. Our extrapolated results were seen to 
agree extremely well with exact results for the spin-half antiferromagnetic XXZ model, up to 
and including the phase transition point at A = 1. Our extrapolated results for the spin-one 
antiferromagnetic XXZ model predicted that the sublattice magnetisation goes to zero at 
A ^ 1.2. This result is in excellent agreement with the value for the onset of the Haldane 



phase of A = 1.167 ± 0.007 ll| for this model. Finally, CCM results for the sublattice 
magnetisation (on both sublattices) of the spin-half /spin-one ferrimagnet were seen to be 
extremely well converged as a function of the truncation index m, and these results were 
finite and non-zero over a wide range of A, up to and including the Heisenberg point at 
A = 1. 

Further applications of the new high-order general-s CCM formalism to the zero- 
temperature properties of lattice quantum spin systems are also envisaged for the future. In 
particular, it is expected that such techniques will be applied to highly frustrated cases with 
spatial dimensionality greater than one which are difficult (if not impossible) to treat using 
other approximate theories. It is also straightforward to extend the ground-state formalism 
)resented here to deal with excited states, as has been done previously for the spin-half case 
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APPENDIX A 

1. Commutation Relations and The High-Order General-s CCM FormaUsm 

In this article we present a new formalism and results for high-order ground-state CCM 
calculations for general spin quantum number, s, based on a model state in which all spins 
on the crystallographic lattice point downwards along the local z-axes. A large part of the 
new formalism relies on the new "high-order" CCM operators defined by Eq. (fT^ and 
also their commutation relations with the single-spin operators in order to determine the 
similarity transforms of various operators, such as the Hamiltonian for example. In order 
to determine these commutation relations we firstly remind ourselves that the ket-state 
correlation operator S is given by Eq. (j21) with C/ = sf^sf^ ' ' ' ^^'^ = Sij^,i2,---,in 
hence 

^ = E E s,,..,-,n«---4 , (Al) 

where each of the indices {ii,i2, ■ ■ - ,11} runs over all lattice sites with the condition that there 
can be no more than 2s of them at any particular lattice site. The usual spin commutation 
relations of the spin operators also apply, 

[s+,s,7] = 2sfV ; [st,s^] = ±sf6i,u . (A2) 

We also note that the commutation of a given operator with S must be distributive, such 
that 

[st,s] = E E '5,,,,.,J[.^,<]4---4 + 4[.^,4]4---< 

+ ••• , (A3) 

where a = {z, +, — }. As pairs of spin-raising operators always commute, we may therefore 
state that [5^,6*] = 0. Furthermore, for the case of [5^,5*] we note again that each index 
runs over all lattice sites, which implies that each term on the right-hand side of Eq. ()A3|) 
is equivalent and that, as there are / such terms, we may write this expression as 

[4, 5] = E E ^ ■ ■ • spt = F,st . (A4) 
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Note again that the "high-order" operators such as Fk are defined by Eq. (fT^. We lastly 
calculate the commutator [s^, S] in Eq. ()A3|) . and using the basic commutation relations of 
Eq. ()A2|) we thus have 

I ii,i2,---,ii 

+ ■■■ , (A5) 



We now commute the operator si past the strings of spin-raising operators in Eq. ()A5|) 
using the basic commutation relations of Eq. ()A2|) . Thus, for example, 

■ ■ ■ 4 = (^k,i2 + Sk,is + ■■■ Sk,i,)sls+ ■ ■ ■ s+ 

+ 44 ■ ■ ■ 44 (A6) 



By inserting Eq. ()A6|) into Eq. ()A5|) we find that 

i-i 

[sk,s] = - 2 ^ iY.)^kAi3,-,ii 44 ■ ■ ■ 44 

~ ^E E ^^k,i2,-,ii 424z ■ ■ ■ "^^"Sfc 

= -Gk,kst - msl , (AT) 

using the definitions in Eq. p3|) . We note again that for the case s = 1/2 the operator 
Gk,k = 0. 

By making use of the nested commutator expansion for the similarity-transformed oper- 
ators [c.f., Eq. (fTUj)]. it is now a simple matter to verify the relations in Eq. (fT^. by using 
Eqs. dXH) and (IXTfl . 

In order to determine the similarity transform of the Hamiltonian however it is also 
necessary to also know the commutation relations of the single-spin operators with respect 
to Fk, F^, Gk,m., and Mk,m,n- The proofs of these commutation relations follow a similar 
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pattern to the proofs given above, and so we merely state them here: 





T 

^k,m^k ' 










^k^ ^ mi 






■Sfc , Fm] 


— '^Gk,mSk 


Mk^k,mSk 5 


^k ' ^rril 


_ 9^2 „+ 
~ '^'^k,m^k 


2FraMk^k,mS'l 




^k,m,n,p'^ k ' 








Nk,k,m.,nSkl 



"k ; 



(A8) 



We note once more that the operators F^, Gk,m, Mk^m,n, and Nk,m,n,p are defined by Eq. 



2. Enumeration of the Fundamental Clusters 



At a given level of approximation, we choose only one of the Nsil^^vi possible symmetry- 
equivalent configurations for a given fundamental configuration of / spin-raising operators, 
where Nb is the total number of Bravais lattice sites, and where vj is a symmetry factor 
dependent purely on the point-group symmetries (and not the translational symmetries) for 
the crystallographic lattice in question and for fundamental configuration /. We note that 
there are Np such fundamental configurations. The first part of the computational algorithm 
is to enumerate all of the "lattice animals" which define the "locale" in which the clusters 
must lie. For the levels of approximation shown in this article it is possible to do this by using 
a simple recursive algorithm which enumerates all possible lattice animals of m contiguous 
sites. This "locale" is explicitly assumed here to be the same for both the LSUBm and 
SUBn-m approximation schemes. Secondly, one then needs to enumerate all possible ways 
in which one can place (2s) or less spin-raising operators on each of the positions of the m 
sites within each of these lattice animals. There are thus (2s)™ possibilities for each lattice 
animal. However, one must also restrict the total number of spin-raising operators to be less 
than or equal to n for the SUBra-m approximation scheme. We note however that there is no 
such restriction on the total number of spin-raising operators for the LSUBm approximation. 
This process thus enumerates all possible connected and disconnected clusters, and we make 
a restriction that we include only those clusters which are inequivalent under the point and 
space group symmetries of both the lattice and the Hamiltonian. A further restriction for 
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the systems under consideration in this article is that we must restrict the set of fundamental 
clusters to include only those which preserve the relationship, = Z^i -^f = 0, with respect 
to the original ( "unrotated" ) Neel model state since [sf., if] = and the ground state lies in 
the = sector. 

3. The Ket-State Equations 

We now wish to determine the CCM ket-state equations, where the J-th such equation 
is given by 

Ei = ^{^CYe~^He^\^) = Q ,V/^0 , (A9) 

where Ai is a normalisation factor given by Ai = (<I>|C7C/|$) = 
• ■ ■ Sj^) {sf_^sf^ ■ ■ ■ sfj\^) . We note once more that we choose only one of 
the NB{U)iyi possible symmetry-equivalent configurations for a given fundamental configu- 
ration in C]~ in order to pattern-match with the terms within H\^) and thus determine the 
J-th CCM ket-state equation. We then computationally match the individual spin-lowering 
operators in C]~ , defined by Eq. (fTH|) . to the spin-raising operators in H\^). We therefore 
put constraints on the indices in the CCM ket-state correlation coefficients, {iSj^^j2^...^j,}, and 
these constraints on the indices allow us to enumerate all possible terms which contribute 
to the CCM ket-state equations. For example, we may consider a specific term in the 
evaluation of the CCM ket-state equations, given by 

{'^\C7 Sirj<l>) = ($1^7 (^sls'^ + F^slsl + Fksts'^ 

+ Gk,mstsi + F,F^stsi^m . (AlO) 

For the case of the linear chain Heisenberg model we let k run over all lattice sites on the 
ID chain and we set m = z + 1. Now consider a specific term within Eq. ()A1Q|) . given 
by ($|C7-Ffc-Fm'5fc We match the indices of the spin-lowering operators in to the 

spin raising-operators in Fk, F^, , and s^. Hence, both k and m are fully constrained 
to take on site values dependent on those indices of the spin-lowering operators in the 
fundamental configuration chosen for C]~ . We note however that one may not always have 
such complete constraints on the indices k and m in H. For example, we may attempt to 
evaluate such a term as, {^\Cj' F^FmSls^l^) . In this case, both k and m are free to cover 
all lattice sites independently from the fundamental cluster utilised in C/. However, we may 
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retain only those configurations in the set {<Sij^,i2, --,ii} in Fk and Fm which are equivalent to 
the fundamental set of configurations under the symmetries of the lattice. This condition 
is sufficient to render the computational problem both tractable and efficient. Finally, the 
resulting coupled, non-linear CCM ket-state equations are easily solved computationally (for 
example, via the Newton-Raphson method) at a given value of the anisotropy parameter A. 



(All) 



4. The Bra-State Equations 

The bra-state coefficients Sj of Eq. (j2)) are formally determined by Eq. (jH)). However, 
this form of the bra-state equations is slightly cumbersome to use, and a simpler and more 
elegant approach is possible by defining the following new set of CCM correlation coefficients 
given by 

xi = Si 

XI = ^SiAMV.) 

where again Ai = ($|C7C/|$) = ($|(s,^s,^ ■ ■ ■ s^^) (s+s+ ■ ■ ■ s+)|$). Note that A^^ is the 
number of Bravais lattice sites. Note also that for a given cluster / then z// is a symmetry 
factor which is dependent purely on the point-group symmetries (and not the translational 
symmetries) of the crystallographic lattice and that / is the number of spin operators. We 
note that the coefficients Ai, z//, and Nb however do not need to be explicitly determined 
because they always cancel out when obtaining ground-state expectation values (see below). 
The CCM bra-state operator may thus be rewritten as 

Np ~ 

S^l+Nj:fCT , (A12) 
1=1 

such that 

H = Eo + Nj^xiEi . (A13) 

We note again that the ground-state energy expectation value is defined by Eq = 
{^e'^He'^l^) and that Ei is the J-th CCM ket-state equation defined by Eq. (jXHI)- The 
CCM ket-state equations are easily rederived by taking the partial derivative of H/N with 
respect to xi, where 

,^Mm^E,. (A14) 

dxi 
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We now take the partial derivative of H/N with respect to xj such that the bra-state 
equations take on a particularly simple form, given by 

= ffl^ = »)+gx.^^ . (A15) 

oxi oxi jr^ OXi 

This equation is easily solved computationally, once the CCM ket-state equations have 
been determined and solved, and the numerical values of the coefficients {xj} may thus be 
obtained. We note that this approach greatly simplifies the task of determining the bra-state 
equations because we never need to explicitly determine the factors A''^, Aj, or z//. 



5. Expectation Values 

Expectation values of spin operators may be treated in an analogous manner to that of the 
expectation value of the Hamiltonian, given by H. For example, the sublattice magnetisation 
for the spin-half and spin-one antiferromagnets of Eq. (|18|) may be written as 

1 ^ ~ 

"5 7=1 i=l 
1 Np 

= 1--J2 m) xixi , (A16) 
^ 1=1 

where i runs over all lattice sites. We again note that the factors Aj and z// in Eq. ()A16|1 have 
cancelled out. Equation ()A16|1 is easily evaluated once the ket- and bra-state equations have 
been solved at a given value of the anisotropy parameter A. (Results for both the ground- 
state energy and the sublattice magnetisation are given in Section III of this article.) 

The situation for the ferrimagnet is slightly different, because the unit cell now contains 
two spins. Thus the magnetisation on A^i the spin-half sites (s = 1/2) is given by, 

1 ^1 
N XT 



sN, f^, Ai 



h=i 



Nf z. Ni 

l-4E7^(*|C7E(^n4)l*) • (A17) 

7=1 ii=l 
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Note that ii runs over all A^i(= N/2) spin-half lattice sites, such that we have the factor 
= 4. It is a simple matter to explicitly enumerate all combinations of orderings of the 
/ spin-raising operators in Z^i^Li -^n which match with spin-lowering operators in , 
although we must also explicitly restrict ii to be a spin-half site. (We note that the factors 
Aj again cancel out with coefficients in {^\CY J2i'i=i{Fhsti)\^) this point.) A similar 
expression may be obtained for M2. 

We may also determine numerical values for other expectation values, such as the spin- 
spin correlation function, given by 

1 ^ ~ 

k=i 

We see from Eq. ()A10|) above that this expression may be written in terms of our high-order 
CCM operators as 

1=1 k=l ^ 

+ Gk,k^rStsUr + FkFk+rStst+)j\^) ■ (Al8) 

The right-hand-side of Eq. ()A18|) may be evaluated computationally in exactly the same 
manner as for the ket-state equations, although no results for spin-spin correlation functions 
are quoted in this article. 
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FIG. 1: CCM results for the ground-state energy per spin of the spin-half XXZ model on the 
linear chain using the LSUBm approximation scheme with m = {2,4,12}. The LSUBm results 
for m = {6, 8, 10, 12} are extrapolated in the limit m ^ oo for this case and are compared to 
exact results of the Bethe Ansatz 1, ^ ^ 4|. (Note that the LSUBm and SUBm-m approximation 
schemes are equivalent for this model.) 
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FIG. 2: COM results for the ground-state energy per spin of the spin-one XXZ model on the linear 
chain using the SUBm-m approximation scheme with m = {2, 6, 8} and the LSUBm approximation 
scheme with m = {2,6}. LSUBm and SUBm-m results are extrapolated in the limit m — >^ cx) for 
A > 1. Note that boxes indicate the CCM critical points. 
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FIG. 3: CCM results for the sublattice magnetisation of the spin-half XXZ model on the linear 
chain using the LSUB?ti approximation scheme with m = {2,8,12}. The LSUBm results for 
m = {6^jl0, 12) are extrapolated in the limit m — > oo and are compared to exact results of Bethe 
Ansatz[l|, y, y, Note that the extrapolated CCM results follow the qualitative behaviour of the 
exact results quite closely and that the extrapolated results go to zero at A = 1.07. 
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FIG. 4: CCM results for the sublattice magnetisation of the spin-one XXZ model on the linear 
chain using the LSUBm approximation scheme with m = {2, 6} and SUBm-m approximation 
scheme with m = {2,8}. Results for the LSUBm approximation scheme with m = {2,4,6} and 
SUBm-m approximation scheme with m = (2, 8} are extrapolated in the limit m ^ oo. Note that 
the extrapolated CCM results go to zero at A « 1.19 and A « 1.18 for the SUBm-m and LSUBm 
approximation schemes, respectively. 
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FIG. 5: CCM results for the sublattice magnetisation of the spin-half spins, Mi, and of the spin-one 
spins, M2, for the spin-half/spin-one XXZ ferrimagnet on the linear chain using SUBm-m approx- 
imation scheme with m = {2,4,6,8,10}. Note that the LSUBm and SUBm-m approximation 
schemes are equivalent for this model at the levels of truncation index m shown. 
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TABLE I: Results obtained for the spin-half Heisenberg model on the linear chain using the CCM 
LSUBm approximation scheme with m = {2, 4, 6, 8, 10, 12}. Np denotes the number of fundamental 
configurations for the ground state. The ground-state energy per spin, Eg/N, and the sublattice 
magnetisation, M are shown. The LSUBm results for m = |6, 8, 10, 12| are extrapo ated in 
the limit m — > cxd for this case and are compared to exact results of the Bethe Ansatz|l|, |2, y, 
^. (Numbers in brackets for the extrapolated CCM results indicate the estimated error in the 
last significant figure shown. Note that the LSUBm and SUBm-m approximation schemes are 
equivalent for this model.) 





Np 


Eg/N 


M 


LSUB2 


1 


-0.416667 


0.666667 


LSUB4 


3 


-0.436270 


0.496776 


LSUB6 


9 


-0.440024 


0.415771 


LSUB8 


26 


-0.441366 


0.365943 


LSUBIO 


81 


-0.441995 


0.331249 


LSUB12 


267 


-0.442340 


0.305254 


Extrapolated CCM 




-0.44315084(6) 


-0.01876(3) 


Bethe Ansatz 




-0.443147 


0.0 
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TABLE II: Results obtained for the spin-one Heisenberg model on the hnear chain using the 
LSUBm approximation scheme with m = {2, 4, 6} and SUBm-m approximation scheme with m = 
{2,4,6,8,10}. Np denotes the number of fundamental configurations for the ground state. The 
ground-state energy per spin, Eg/N, the sublattice magnetisation, M, and the critical values of 
Ac are shown. The SUBm-m results for m = {4,6,8} and LSUBm results for m = {2,4,6} are 
extrapolated in the limit m ^ oo in this case, and are compared to results^^ of the DMRG method 
P]. The ** symbol indicates that extrapolated CCM results for the sublattice magnetisation, M, 
go to zero at = 1.19 and A^ = 1.18 for the SUBm-m and LSUBm approximations, respectively, 
also given in this table) and are thus ill-defined below these points. The result* for the position 



of the phase transition point using the large-cluster-decomposition Monte Carlo method 
quoted. 



ll| is also 





Np 


Eg/N 


M 


Ac 


SUB2-2 


1 


-1.316625 


0.809068 




SUB4-4 


7 


-1.360084 


0.694610 




SUB6-6 


37 


-1.375607 


0.607339 


0.593 


SUB8-8 


247 


-1.383466 


0.533252 


0.249 


Extrapolated SUBm-m 




-1.408039 


** 


1.19t 


LSUB2 


2 


-1.320608 


0.800702 




LSUB4 


11 


-1.369428 


0.646536 




LSUB6 


63 


-1.383292 


0.532198 


-0.670 


Extrapolated LSUBm 




-1.403737 


** 


1.18t 


c.f. 




-1.401484038971(4)tt 


0.0 


1.167(7)* 
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TABLE III: Results obtained for the spin-half/spin-one ferrimagnetic HAF on the hnear chain using 
the CCM SUBm-m approximation scheme with m = {2, 4, 6, 8, 10} compared to results of DMRG 
calculations. CCM results for m = {6, 8, 10} are extrapolated in the limit m oo. Np denotes 
the number of fundamental configurations for the ground state. The ground-state energy per spin, 
Eg/N, the sublattice magnetisation on the spin-half sites, Mi, the sublattice magnetisation on the 
spin-one sites, M2 are shown. (Note that the LSUBm and SUBm-m approximation schemes are 
equivalent at the levels of approximation shown for this model.) 





Np 


Eg/N 


Ml 


M2 


SUB2-2 


1 


-0.70710678 


0.70710678 


0.85355339 


SUB4-4 


5 


-0.72582592 


0.59865621 


0.79932811 


SUB6-6 


21 


-0.72697237 


0.58611255 


0.79305628 


SUB8-8 


93 


-0.72704344 


0.58503667 


0.79251834 


SUBlO-10 


427 


-0.72704696 


0.58497261 


0.79248630 


Extrapolated CCM 




-0.7270474 


0.5849641 


0.7924820 


DMRG 




-0.72704 


0.58496 


0.79248 
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